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Abstract 

This is the first part of a threefold article, aimed at solving numerically the Poisson problem 
in three-dimensional prismatic or axisymmetric domains. In this first part, the Fourier Singular 
Complement Method is introduced and analysed, in prismatic domains. In the second part, 
the FSCM is studied in axisymmetric domains with conical vertices, whereas, in the third part, 
implementation issues, numerical tests and comparisons with other methods are carried out. 
The method is based on a Fourier expansion in the direction parallel to the reentrant edges of 
the domain, and on an improved variant of the Singular Complement Method in the 2D section 
perpendicular to those edges. Neither refinements near the reentrant edges of the domain nor 
cut-off functions are required in the computations to achieve an optimal convergence order in 
terms of the mesh size and the number of Fourier modes used. 
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1 Introduction 

The Singular Complement Method (SCM) was originally introduced by Assous et al 
(HI El) for the 2D static or instationary Maxwell equations without charges. The cases 
with charges have been recently solved by Garcia et al jHlE]) including the numerical 
solution to the 2D Vlasov-Maxwell system of equations. The SCM has been extended 
in to the 2D Poisson problem. Further extensions to the 2D heat or wave equations, 
or to similar problems with piecewise constant coefficients, can be obtained easily. As 
a matter of fact, this stems from the analysis which is performed hereafter (see Re- 
mark The primary basis of the SCM is the decomposition of the solution into 
regular and singular parts. Methodologically speaking, the SCM consists in adding 
some singular test functions to the usual Pi Lagrange FEM so that one recovers the 
optimal //^-convergence rate, even in non-convex domains. In 2D, one may simply add 
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one singular test function per reentrant corner. 

There exist a couple of numerical methods in the literature for accurately solving 
2D Poisson problems in non-convex domains. It was shown in [T^j that the SCM can 
be reformulated so that it coincides with the approach of Moussaoui I 27( ^ when L- 
shaped domains are considered. The SCM differs from the Dual Singular Function 
Method (DSFM) of Blum and Dobrowolski JHI in that it requires no cut-off functions. 
Actually, when the numerical implementation of the SCM is carried out, the cut-off 
function is traded for a non-homogeneous boundary condition. Note that Cai and Kim 
|12j recently proposed a new SFM which involves the evaluations of singular and cut- 
off functions and the solution of a nonsymmetric elliptic problem. The SCM is clearly 
different from (anisotropic) mesh refinement techniques |28[ IH} l2^ 1^ [2]. and can be 
applied efficiently to instationary problems (see Remark 14. 1|) , since it does not need 
mesh refinement and thus larger timesteps may be allowed. However the anisotropic 
mesh refinement methods have one advantage: they require only a partial knowledge 
of the most singular part of the solution. 

The numerical solution of 3D singular Poisson problems is quite different from the 
2D case, and much more difficult. This is a relatively new field of research: most 
approaches rely on anisotropic mesh refinement, see for instance pi 151 OH HI I25j. and 
|2j and Refs. therein. To our knowledge, this series of papers is the first attempt to 
generalize the SCM for three-dimensional singular Poisson problems. Specifically, we 
shall consider the numerical solution of the Poisson problem: 
Find u G HQiO.) such that 

-Au = f in n, (1) 
where / S L^(0), and Q is a (right) prismatic domain described by 

Vt = uJxZ, (2) 

and (J is a two-dimensional general polygonal domain, Z is an interval varying from 
to a positive constant L on the xs-axis. The bases of the domain are the subsets of the 
boundary dVt^ which are included in the planes {x3 = 0} and {x3 = L}. 

The case of an axisymmetric domain is considered in the companion paper jl4) . 
When the Poisson problem is solved in this class of domains, two difficulties arise. 
The first difficulty is that one has to deal with weighted Sobolev spaces, the weights 
being functions of the distance to the axis. The second one is that there exist two 
kinds of geometrical singularities: reentrant edges like in the prismatic case, and, in 
addition, sharp conical vertices. As for implementation issues and comparisons with 
other methods (such as variants of our method, the FSCM, or mesh refinement tech- 
niques ^1), we refer the reader to jl5j . 

The rest of the paper is organized as follows. In the next Section, some theoretical 
results concerning the regularity of the solution to the Poisson problem in prismatic 
domains are recalled. A priori regularity results of the solution u to and a first 
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splitting of the solution into regular and singular parts, are emphasized. In Section 3, 
some results about the Fourier expansion along X3 are recalled and/or proven. This 
suggests a framework for building the Fourier Singular Complement Method (FSCM) 
for accurately solving the problem ^ , using a Fourier expansion in 2:3 , and an improved 
variant of the Singular Complement Method in the 2D section uj. In Section [IJ 
we study the variant of the SCM, based on a theoretical splitting of the solution 
to 2D problems of the form —Au^ + /xu^ = in u; (with a parameter /i > related 
to the Fourier modes). The main feature of the regular-singular splitting is that it 
is chosen independently of /i; this independence is important, and very helpful, from 
the computational point of view. Estimates on Sobolev norms of and its splitting 
are established. To end this Section, the SCM is considered from a numerical point 
of view, to approximate accurately, via the discretization of the splitting: the opti- 
mal H^-norm convergence of the order 0{h) is recovered. In the last Section, we first 
prove a refined splitting of the solution u to the 3D Poisson problem under suitable 
assumptions on the right-hand side /, using the Fourier expansion along X3. Then, we 
build the numerical algorithms which define the FSCM, and we show that the FSCM 
has the optimal convergence of order 0{h + N~^), where h is the 2D mesh size and N 
is the number of Fourier modes used. 

Throughout this paper, when two quantities a and b are such that a < C b, with a 
constant C > which depends only on the geometry of the domain, we shall use the 
notation a < 6. 

2 Poisson problem in prismatic domains 

Let us recall that a (right, open) cylinder of M^, with axis parallel to X3, is equal to 
Dx I, where D is any connected (open) subset of M?, and I is any (open) interval of M. 
Let us proceed then with some remarks on the class of domains 0, i.e., the prismatic 
domains. A priori, such domains could be considered: 

• either as truncated infinite cylinders; 

• or as polyhedra. 

As it happens, considering Q as a polyhedron is helpful, in a simple manner. Indeed, 
from IHEI) we know that, in any polyhedra, the solution u to Q can be split as 

U = Ur + Ue + Uy, with Ur G -ff^(O), (3) 



Above, Ur is called the regular part, Ug the edge singularity part, and Uy the vertex 
singularity part. Note that when Ue ^ or Uy ^ 0, they do not belong to H^{i^). 
The summation in Ue is taken over all reentrant edges e, (pe, (/"e, -^e) denote the local 
cylindrical coordinates, and Tr/ue the dihedral angle (so that ae E]1/2,1[). Last, the 
summation in Uy is taken over all non- convex vertices v and over all eigenvalues A„ of 
the Laplace-Beltrami operator, which belong to the interval ] — l/2,l/2[, and (p^, 6y,(j)y) 




l/2<A„<l/2 





3 



denote the local spherical coordinates. 

In our case, i.e., when $7 is a prismatic domain with polygonal bases, it has been shown 
m that the vertex singularity part always vanishes, so (jSl reduces to 

u = Ur + Us, with Ur £ H'^{i^) and Us = ^J^eiPe, Ze) sin(ae0e). (4) 

e 

Let us describe how one can fall into the other class, that of the infinite cylinders. 
The first step is to introduce a suitable continuation u of the solution u (odd reflection 
at the bases) along the X3 direction from Z to M: one builds a problem to be solved in 
the infinite cylinder Coo = w x M. Unfortunately, with this continuation technique, 
one gets a solution (and data) which is not in L^(Coo). Thus, one introduces in 
a second step a smooth truncation function rj, such that rj{xz) is equal to one for 
X3 s] — i/2, 3L/2[, and to zero for |x3| > 2L. Then, one multiplies u by r/, to obtain a 
Poisson problem in Coo with solution = wq. This time, one has G H^{Coo) (and 

= —Am'' G L^(Coo)). By construction, the restriction of on J7 coincides with u. 

Interestingly, it has been proven in |211l26j . that a splitting similar to (jlj holds for 
u'^. Furthermore, Ms can be expressed as 

Us = leiPe, X3)Pe" sin(ae(/)e). (5) 

The function 7e in ((5} is often called in mechanics the stress intensity distribution. 
On the one hand, in the original paper 7e is expressed as a convolution product. 
On the other hand, in it is characterized as the solution to a second order PDE. 
Finally, the regularity of the singular part m^, can be expressed accurately as follows 
|26| . Let 6 denote the minimal distance between two reentrant edges, and for each 
reentrant edge e, let ^l^ = {x £ Q : d{x,e) < 6/2}. Then 

' M G H^+''~%n)_, Ve > 0, a = miUeOe, 

uGH''{n\Uene), 

ptdiu E H^i^e). Ve, V/3e > 1 - ae, ^ = 1, 2, ^ ' 

^ dsuGH^n). 

In Sectional using a Fourier expansion along the X3 axis, we recover some proper- 
ties which are very similar to (^HHl). 

We end this Section with remarks on other possible boundary conditions. 

If the boundary condition for u on the bases of the physical domain 0, are the 
non-homogeneous Dirichlet boundary condition: 

u = g at 2:3 = and X3 = L, 

one can set w = u — g with g being a continuation of g into fi. Then the problem 
reduces to the case with the solution w satisfying the homogeneous Dirichlet boundary 
condition, assuming that g e H^{n). 
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If the boundary condition for u is the homogeneous Neumann boundary condition 
dnU = on f?i7, then one can replace the sin(ae0e) factor in Q by the expected 
cos(ae(/)e). Moreover, to obtain an expression hke one uses an even reflection of 
u at the bases of the domain. If we have the non-homogeneous Neumann boundary 
condition: 

dnU = g at 2:3 = and X3 = L, 
one may then study the solution 

w{x) = u{x) — / g{xi,X2,z)dz 
Jo 

first, which satisfies the homogeneous Neumann boundary conditions at the bases of 
the domain 0. Here ^ is a continuation of g into which is again assumed to belong 
to H^{n). 

From now on, we assume, for ease of exposition, that the polygon lo has only one 
reentrant corner C, i.e., with an interior angle larger than vr, denoted as vr/a, with 
1/2 < a < 1. In particular, the summation which defines the singular part Us in 
reduces to exactly one term. 



3 Fourier expansion 



We devote this Section to some justifications about the Fourier series expansion of the 
Poisson solution to First, one can show, following for instance Heinrich's proof of 
Lemma 3.2 in the well-known result 



Lemma 3.1 For any f G there exist Fourier coefficients defined by 

2 k-K 
fk{xi,X2) = -j /(xi,a;2,X3)sin— xsdxs, /c = 1, 2, 3, • • • , (7) 



such that fk £ L {ijo) and 



00 ^ 
/(a;i,X2,X3) = ^/fe(a;i,X2)sin— X3 a.e.inQ, (8) 



and 



L 

k=i 



L °° 



2 

L2(f7) — ^ 

k=l 



k=l 

If f e Hl{Vt), then fk G Hq{uj) for all k and 

iiv/iii.(n) = 2 E {ii^/'^iii^H + (t) ii^'^ii'^h} < (10) 

k=l 
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For / in L'^{Q), let us introduce the sequence of partial sums {Fk)k of the Fourier 
decomposition of /, which converges to / in L'^{Q), cf. ©: 



K 



FK = Y.fk 



k=l 



sin — x^, for K > 0. 
L 



We note that when / is in Hq{Q), {Fk)k converges to / in Hq{Q), according to (jTUIl. 
Also, the sine functions can be replaced by cosine functions with the same argument 

^X3, and (Tfrmil still holds (for with any / in H^{VL).) 



L 



In our subsequent analysis, summations like 



k\\L'2{u,) 



(12) 



k=l 



will appear. The result below provides a characterization of elements / of L^(r2), which 
are such that (fT^ is bounded. Let the following Sobolev spaces be introduced: 



hUn) 



H\]^,L[,L\u))={f eL\n) : ^3feL\^l)}■ 
H',{]0,L[,L\u;))={feh\n) : /||,3=o} = /|{.3=L} 
H\]0,L[,L\u;))={f Gh\n) : d^sf G L^Q)}. 



0} ; 



Lemma 3.2 Given f G L'^{Q), one has the following equivalences 

oo 

/e/ii(l^)^^A;2||/,||2,^^^<oo; (13) 

k=l 

oo 

/ G hliQ) n h\n) ^Y^k^hWl,^^^ < oo. (14) 

k=l 

Proof. Let / be in L^{n). 

Assume in addition that / G hl{Q). We note that, by the definition of the Fourier 
mode fk and integration by parts (/ vanishes at the bases), one has 

{kTr)fk = ~2 y / "^^3 = 2 J Sa/cos— xadxs. 

oo 

Since by assumption, ^3/ is in L'^{Q), one gets the expected ^^Il/A:|li2((^) < o<d. 

k=l 

Let us prove the reciprocal assertion. For / in L^(J7), as the sequence {Fk)k (see 
(fTT|) ) converges to / in L^(r2), one infers that (^s-Fa')^' converges to ^3/ in H~^{^). 
Now, if the sum is bounded, {d^FK)^ is a Cauchy sequence in ^^(0), so it converges 
in this space, and its limit ^3/ is in L^(Q). Since both {Fk)k and {d-iFK)K converge 
in L'^{Q), {Fk)k converges to / in h^{Q), and, as Fk belongs to h\{Q) for all K, f is 
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also in hK^l), which proves (fT^ . 

In order to establish H14|). one proceeds similarly, by performing a second integration 
by parts. Note that for this additional integration by parts, no assumption is required 

on the trace of / at the bases, since the (sin —x^)k vanish there. 

Jk = ^-r osf cos— X3dx3 = -2 dsa/ sm— xadxs. 



L L Jq L ./n L 

With this identity, one concludes the proof easily. 



Now that the general results have been obtained, we focus on the Poisson prob- 
lem . Consider the weak form of the Poisson problem: 

a{u,v) = f{v) \/veH^{n) (15) 

where a(-, •) and /(•) are given by 

a{u,v) = / Vu-Vvdx, f{v)= / fvdx. 
Jn Jn 

We expand the solution u in in the Fourier sine series: 

oo 

n(xi, X2, xs) = ^ Ukixi,X2) sin -j^xs. (16) 

k=l 

Following again Heinrich's proof of Lemma 3.2 in j23j . the next two Lemmas hold. 
Lemma 3.3 For any u,v £ HQ{i^), we have 

^ oo Ij 

a{u, v) = - ak{uk,Vk), f{v) = ^ X] fk{vk), 

k=l k=l 

where and fk are given by 

ak{uk,Vk) = j jVufc • Vvk + UkVk^dxidx2, fk{vk) = j fkVkdxidx2, 

and Uk, Vk and fk are Fourier coefficients of u, v £ i/o(r2) and f G L'^{Q) respectively. 

Lemma 3.4 For any f G L'^{Q), let u G Hq{Q) be the unique weak solution of US]) 
and Uk and fk be the Fourier coefficients of u and f. Then Uk G Hq{lo) is the unique 
solution of the following 2D weak problem: 
Find Uk G Hq{lo) such that 

ak{uk,v) = fkiv) yvGHQitu). (17) 
Moreover, Uk satisfies the following a priori estimates: 



l^{\^Uk\^+{'^yul]dx,dX2 < 



A; = 1,2,--- 



k=l 
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This means that the fc-th Fourier mode of u is characterized as the unique solution to 

the 2D problem 

Find Uk G Hq{uj) such that 



-Auk + ly'Y' j = fk in Lo; Uk = on dto. (18) 

As Corollaries, one gets a convergence result of the sequence of partial sums {Uk)k 
of the Fourier decomposition of u, and also the last result of ©• 

Corollary 3.1 Let f G L^(r2), and u be the solution to Then {Uk)k converges 
to u in H^{il), and {AUk)k converges to —f in L^(r2). 

Proof. The fact that {Uk)k converges to u in H^{Q) is a consequence of Lemma l3. II 
Then, one notes that 

-AUk = 2_^{-Auk+ {— j Uk) sm—X3 = 2^ fk sm—xs = Fk, 

k=l ^ ^ k=l 

which yields the result on the convergence of {AUk)k- o 



Corollary 3.2 Let f G L'^{Q), and u be the solution to (0^. Then d^u G H^{Q). 
Proof We prove that, for i = 1, 2, 3, di^u belongs to L^(J7). 

oo 

For i = 3, thanks to the last bound of the Lemma [3.41 there holds ^^II^A:|li2(j^) < oo. 

k=l 

Result (O yields u G hl{n) n h'^{n), so that ^331* is in L'^I^Q). 
For i = 1,2, we note that 

disUx = ^ kdiUk cos ^^3. 

k=l 

00 

According again to the last estimate in Lemma 13.41 ^^||9j^A;||/,2(^) < 00, so di^u is 

fc=i 

inL'^in). o 

To conclude this Section, we note that the Fourier expansion ()16(1 of u together with 
the series of 2D problems (jlSp suggest the numerical approximation scheme below, i.e., 
define the Fourier SCM (FSCM) approximation of the solution u to (fT3)) as follows: 

U^{xi,X2,X3) = ^ u^(xi, X2) sin (19) 
fc=i 

where N is the total number of Fourier modes used in the approximation, and is a 
suitable approximation of u^, to be studied in the next two Sections. 
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4 Regular- singular decomposition in the 2D do- 
main uj: theoretical study 

The main interest of this paper is to propose some efficient numerical method for solving 
the three-dimensional singular Poisson problem in a prismatic domain. Basically, 
the method reduces the 3D problem into a series of 2D Poisson- like problems, see H18() . 
by the Fourier expansion of the 3D solution along the xs-direction. This Section will 
thus focus on the 2D singular Poisson problem: 
Find G Hq{lo) such that 

-Au^ +^n^ =/ in (20) 

In the case of the Fourier expansion, one considers /i = k'^ir^/L'^ and / = in 1)21) |). 
Due to the presence of the Fourier mode index k, the coefficient fi varies in a large 
range, from -k'^/L'^ to N'^tt'^ /L?, where is the number of Fourier modes required 
subsequently in the numerical approximation (cf. Section EJ. This brings in one of the 
main difficulties in the subsequent error estimates, which should hold for all ^'s in a 
large range. 

As a preliminary remark, we note that, according to the most singular part of 
the solution to ()2U() is of the form /)" sin(a^), compared to Q) in 3D. 

Let 7i, 72, • • • , be the line segments of (9a;, where 71 and 72 are two line segments 
which form the single re-entrant corner of to. Our numerical method is based on the 
following important decomposition of the space L'^{io) |22j : 

L\uj) = A[H\uj) n H^{uj)] e N , (21) 
where A'^ is a space of singular harmonic functions defined by 

N={p€ L\u) : Ap = 0, p\^, = in {Hll,\^k))\ 1 < A; < i^}. 

Above, the space ifQQ^(7fc) is made up of elements of H^^'^{'yk), such that their continu- 
ation to du) by zero belongs to H^^'^{duj). Its dual space is denoted by {Hq(^ {^i^))' . To 
understand that the boundary condition on p holds in this dual space, let us mention 
that one can prove that, given any 4> in H'^i^) ^ ^oi^)^ 9n4>\'yf. belongs to -frgQ^(7fc). 
Then, the fact that p\^^ = simply reflects a surjectivity property, which states that 
the mapping (f) 1— > is onto, from H'^{uj) Pi Hq{uj) to HqI^"^ (jf.). 

As the domain co has only one re-entrant corner, we know dim(A^) = 1, and 
N=spaji{ps} for some ps £ N \ {0}, see Grisvard [Uj. 

Let (ps be an element in Hq^lo), which solves the Poisson problem 

— A(j)s=ps in oj . (22) 
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Then by the decomposition H21|) . we can spht the solution to equation H2U|) as 

Ufj. = Ufj. + c^(l)s, (23) 

where £ H'^il^) ^ Hq{uj)^ and is called the regular part of u^. 

We will devote the rest of this Section to the derivation of some a priori estimates 
for the solution n^, its regular part and the singularity coefficient c^, as well as the 
solvability of and c^. Let us first introduce some notation. 

Throughout the rest of the paper, uq will be a frequently used fixed positive constant 
lying in the interval ]5,a[, where a s]^, 1[ is the singularity exponent. | • |s is used to 
denote the semi-norm of the Sobolev space H^[lo) for any s > 0, (•, •) and || • ||o are 
used to denote the inner product and the norm in the space L'^{uj). Also, (•, •) will be 
used for the dual pairing between the space -ffo(w) and H~^{uj) when necessary. 

The following lemma summarizes some a priori estimates on and c^. 

Lemma 4.1 Let he the solution to the Poisson problem \2U\) . then we have the 
following a priori estimates: 

/^II^^mIIo < ll/llo, VJ^K\i < ^ll/llo, IIAn^llo <2||/||o, (24) 

|cmI < M^'^ll/llo (25) 
K\i+ao < /i-^ll/llo. (26) 
Proof. Multiplying equation 1)20(1 by and integrating over oj yield 

Wt^ll + M ll^^^tllo < ll/llo \Wf,\\o , 

which proves the first estimate in (|24|) . Then applying the Cauchy-Schwarz inequality, 
we further obtain ^ ^ 

\u^,\l + M II^^mIIo < ll^/^llo + ^ll'^llo' 

which leads to the semi- norm estimate in ()24() . 

The last estimate in ()24|) follows immediately from Au^ = ^u^ — f and the first 
inequality in 

As far as ()25|) is concerned, it is a simple matter to check that the singularity 
coefficient c^, multiplied by some constant /?*, equals the singularity coefficient c(//) of 
jm pp. 62-69]. Indeed, in Grisvard's papers, is decomposed into: 

7x^ = n^ + c(^)e-v^''e(p)p"sin(a</.), e H^{i^) D H^{u,) (27) 

where ^ is a smooth cut-off function, equal to one in a neighborhood of 0. 

On the other hand one can decompose the singular part in (|23() as (cf. or ()46() 
below) 

C;.</.s = c^(<A + /5V"sin(a(^)) , ^ e H\io), (3^ = ^\\p,\\l 
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Using this, (^5)) and ((?7|) . we can write 

(c^r - c(/i)e(p))p° sin(a<^) 
= Ufj,- (n^ + Cfj.(t>) - c{ii)i{p)p°' sin(a(/)) 

= u'; + c{p)(e-^P-l)i{p)p^Bin{a4>)-{u^ + c^^). (28) 

Noting that each term on the right-hand side of (|28() belongs to H'^{lu), we must have 
C/i = c{ij)/(5* . But it is shown in ineq. (2.5.5)] that 

\c{p)\<p-'^\\fh. (29) 

which imphes (|25jl . 

In order to derive the estimate , we shall use (|27i;-i()j) , with the additional norm 
estimate |221 ineq. (2.5.4)] on the regular part u^, namely 

\n';\2 + ^\u%+p\\u%<\\f\\o. (30) 
Indeed, from the estimates 

k^|i<^-'/'ll/llo, |n^|2<||/||o, 
we have then by standard interpolation theory that 

Iw^ll+ao <^"^ll/l|0. 

Next, we use (1^ and a direct estimate of the i/^+^o 

semi-norm to bound the singular 

part in (P7|). Actually, there holds 

J X^UJ J x'duj ^ I 

Due to the uniform smoothness (in p) of e~v^''^(p)/j" sin^acj)) for p > po > 0, it is 
possible to evaluate the integrals only on = {{p-,(t^) S]0, /Oo[x]0, 7r/a[}. Then, one 
performs the changes of variables s = y/pp, s' = ^/pp' , to find 

|e-v^^^(p)p"sin(a(/>)|^i+.o(.^) < C{ao)p-^. 

This with ^ leads to o 
Now, let us study the solvability of and in decomposition H23|) . For conve- 
nience, we introduce the notation a^(-, •) and the norm || • ||a: 

af^{w, v) = {Vw, Vv) + p {w, v) , \\v\\l = a^{v, v) , 

and the linear mapping from Hq{uj) to H^^{uj), defined by A^u = —An -|- pu, or 
equivalently by 

H-^c,) < \w,v >//iH= o.^^{w,v) \fw,v e H^{uj). 
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It is not difficult to verify that is a one-to-one and onto mapping, so it is invertible. 
So, we claim that and solve the following coupled system: 

a^{u^,v) + c^a^{(j)s,v) = {f,v) \/v £ Hq{u;) , (31) 
{WPsWo + ^I'/'sIl) + M {Uf„Ps) = {f,Ps) ■ (32) 

In fact, by multiplying the equation H2U|) by ps and integrating over uj we obtain 

-(Au^,Ps) + {UfM,Ps) = {f,Ps) , 

then ()32() follows readily from the decomposition ()23() . the orthogonality between ps 
and Au^, along with the relation (|22)) and its following direct consequence 

|,/.,|2 = (,/,„p,). (33) 

On the other hand, the solution of H20|) also satisfies the weak form: 

(Vu^,Vf) + /i(u^,?;) = {f,v) yveH^iio). 

This and the decomposition (pS)) lead to the equation (|^T|). 
Below, we show the well-posedness of the system (|3T|) - (|32|) . 

Lemma 4.2 There exists a unique solution (u^^c^) to the coupled system 
and the following stability estimates hold: 



\cu\ < 2 - — — , 2 < 4 

WPsWo 

where Cp is the constant in the Poincare inequality. 

Proof. To see the unique existence, we rewrite ()31() as the following operator form: 

A^Uf, + A^cPs = / in H~\uj). (34) 

As the inverse of A^ exists, we know from ()34() that can be determined if is 
available: 

^1^ = ^Ji^f - Ci,4>s- (35) 
This is exactly our original decomposition (|23|). Substituting this into (|32|) . 

(IIPsIIo + /^l</'s|?)c^ + ;^(^;;V - Cf,(j)s,Ps) = if,Ps) ■ 
With (jnSl), we obtain that 

if-fiA-^f, Ps) 

= ^^^2 • 
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With uniquely determined, is clearly uniquely determined by (|^H) or ()35|1. 

Next, we derive the stability estimates in Lemma [4.2l We show that these estimates 
are the consequences of ()35ll36p and of the following inequality 

\\A-^9h<-\\9h VffGL2(c^). (37) 
A* 

In fact, if 1)37(1 is true, then the desired estimate on follows from (|36|) : 

+ /^PmVIIo 



\cu\ < — — — < 2- 

WpsWo WPsWq 

On the other hand, we have from (|33|) and the Poincare inequality that 

UsWo < Cp\\V<Ps\\o < Cj,\\ps\\o. 

Using this and the bound of c^, we derive from ()31j) by taking v = that 

llVn^llo + ^ll-U/xllo < ll/l|o||ttMl|o + \c^i\ (llVc/fsllollVit^llo + /i ||0s||o||ttMllo) 
< ||/||o||u;.||o + 2Cp||/||o||V7i^||o + 2f,Cj.\\f\\o\\uJo ■ 

Then, an application of the Young inequality yields 



• 



This implies 

l\\nX < Q + 2C|, + 4mC^) ||/||g< (^-L + 2y;[IC^y||/||2, 

so the desired estimate on ||n^||a follows. 

We now show the ff^-norm estimate. By the decomposition (|35|) . we have = 
Aji^f = Uf, + c^4>s, and 

-Aufj, = -Au^ + c^Acps = f - iiu^- c^ps, 

which gives 

||Aii^||o < ll/llo +/i||u^||o + |c^| \\Psh- 

But we know from Lemma |4.1l that /x||u^||o < ||/||o- This, along with the previous 
bound for c^, leads to 

IIAu^llo <4||/||o. 

Now, for any v G H^{io)'^ such that • r = on dio, with r the vector tangential to 
duj, it is well-known (cf. ^2j) that (since a; is a polygon) 



Yl W^kViWl = \\cmlv\\l + ||divi;||g. 



l<fc,K2 
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So, by taking v = Vu^, one actually finds 

\Uf,\2 = IIAn^llo < 4| 

Finally, it remains to prove (|37() . By the definition of a^(-,-), we easily see the 
following lower bound: 

H-^cv) < \v,v >Hi(^^^= a^{v,v) > fi\\v\\l \/veH^{uj). (38) 

Then for any g £ L'^{uj) C H~'^{lo), let v = g G HI{ijo). One has A^v = g in L'^{uj) 
and it follows from H38|l that 

P;'5llo = Ikllo < -(^M^'^) = \{9.A-^g) = i||5||o||A;iff||o, 
ft fi fi 

which proves (|57|) . o 

We end this Section with a number of important remarks on the theoretical and 
practical range of the splitting into regular and singular parts. 

Remark 4.1 Equation \2U\) is also useful when the 2D heat equation is considered in 
uj: 

du 

— -Au = fin ujx]0,T[, 
dt ^ ^ 

with initial condition and (homogeneous) Dirichlet boundary condition. As a matter 

of fact, assume it is first discretized in time, with a time-step 6t, at times tm = m5t, 

m = 0, 1, • • • ; let u™ = u{tm). Then one has to solve in space the implicit problems 

(with 9 g]0, 1] given) 

Find 14""+^ e -H'd(t^) such that 

-An"+i + -rV'"™^^ = f(t +i1 + ^-^f(t ) + tV-"™ + ^^Au™, in uj. 

Above, 9 = 1 (resp. 9 = 1/2) corresponds to the implicit Euler (resp. Crank-Nicolson) 
scheme. This is precisely h20\) with // = l/95t. 
Clearly, implicit schemes for the 2D wave equation 

— - Au = f in Lux]0,Tl 
lead to other instances of Equation \2U\) . 

Remark 4.2 Both (ps and ps in i22\) are chosen independent of fi, f and u^, so their 
norms will be regarded as some generic constants (i.e., independent of fL, f and u^.) 

Remark 4.3 Instead of the decomposition \2S\) . it seems more natural \21\ \2 2^ to take 
the decomposition = u'^ + c^cj)^, where (j)^ S -f^d(w) depends on the parameter fi, and 
it is the solution to the problem: — Ac/)^ + = ^^^^ Pt^ ^ ^fj. \ {0}; where 

is given by 



[p£L\uj) : (-A + ^I)p = 0, p\^, = in (Hli^jk))' , 1 < k < kY 



But the decomposition i2cl\) has an important advantage: the singular part (ps is in- 
dependent of the parameter /x. As we shall see, this will be much less expensive than 
using the above more natural decomposition. 
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5 Discrete formulation in the 2D domain w. the 
SCM 



In this Section we shall formulate the generalized SCM for solving the coupled sys- 
tem (|3i p - (|32p and derive the error estimates of the approximate solutions. The SCM 
was first introduced by Assous et al |8j for solving the 2D static or unsteady Maxwell 
equations without charges, and then used in for the 2D Poisson problem. As we 
will see, the formulation of the SCM for the 2D Poisson-like problem H18|) is quite 
different here due to the involvement of the parameter /i. 

Let 7/i be a regular triangulation of the domain w, with vertices {Mj}^^^^ and 
the last Ni, vertices lying on the boundary duj. We define to be the continuous 
piecewise linear finite element space on Th with the standard basis functions {V'j}^i'^' 
(cf. JH])- We further define Vq to be the subspace of with all functions vanishing 
on the boundary of w. The interpolation associated with the space will be denoted 
by n;,. 

5.1 Approximation of the singular function ps 

We start with the finite element approximation of the singular function ps ^ N m (|22)) . 
Recall the splitting (see 

Ps=P+Pp, p£H^{uj), Pp = p~°'sm{a(l)) . 

As Ps is harmonic in uj, one can directly verify that the regular part p in the splitting 
solves the problem: 

Find p £ H^{uj) such that p = s on dtu and 

( Vp, Vv) = V w G (cj) (39) 

where the boundary function s is given by 

s = on 7i U 72; s = —pp on 7a; (3 < < K). 

For the finite element approximation of the problem (|39j) . we shall use the simple 
treatment of the boundary condition: 

7^hS= Yl (40) 

Then we approximate ps by p^ = ph +Pp , where ph is the piecewise linear finite element 
solution to the problem (|39|) . Namely, Ph = T^hS + p^ where p^ G solves 

(Vph,Vvh) = yv^eVo^. (41) 

The error estimates for the singular function ps and its finite element approximation 
Pg are summarized in the following lemma. 
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Lemma 5.1 We have ^ 

\Ps-Ps\i<h"\ \\Ps-p%<h^'''- 
Proof. We introduce a smooth continuation of s into ui: 

s = -p^{l-C{p)). 

Clearly, s = s on du) and s S H'^{lo). Let = p — s. It is known that p S H^~^°'°{io), 
so we have p° e n H^iu). It follows from ^ that 

(V/,Vi;) = -(Vs,V?;) Vt;G/7o^(w). (42) 

Recall Hh is the interpolant associated with V^, thus we can rewrite the finite element 
solution ph to the system (|1T|) as ph = + P^ with p^ e Vg^ now solving 

(VplVvh) = -{VIlhS,Vvh) yvheVo\ (43) 

by noting IlhS = iTh s on dco. 

Now we are ready to derive the error estimates. It is clear from (|42j) and (|4.S|) that 

(V(p°-p^),Vi;;,) = (V(n;,s-s),V^;^) V^;;,Gyo''- (44) 

Using this, we obtain for any G Vq that 

II V(/ - g,)f > II V(/ - p°)f + 2(V(n,5 - 5), V(pO - q,)), 

taking q^ = HhP^ above and using the Young inequality leads to 

\p'-pI\1 < |/-n^p°|? + 2|n,5-5|i(|p°-/|i + |/-n;,/|i) 



< 2 1/ - Uhp^'ll + n\Ph- p°\l + 3 - ~s\i . 



Then by the standard interpolation results we obtain 

\p' - P^l? < 4 |p° - u,pX + 6 |n,.~ -s\l< /.2-|p°|?+„„ + \rs\l 

This leads to the desired ff^-norm error estimate: 

\Ps-Ps\l = \P - Ph\l = \P° + S - pI - IlhS\i 

< 1/ - P°|l + N~ - n;.s|i < + h\s\2 < 

Finally, we apply the Nitsche trick to derive the L^-norm error estimate. Let 
w £ Hq{uj) be the solution to the variational problem 

{Vw,Vv) = {p'' -pIv) yveH'oiu;). (45) 



^By construction, neither ps nor p'^ belong to H^{u!), due to the presence of Pp, but the following holds: 
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By the elliptic theory, we know w £ 11^^°^° (uj) and 

Hl+ao < llf° -P°l|0- 

Let Wh be the finite element approximation oi w: Wh £ Vq solves 

{vwh, vvh) = - plvh) yvhe Vo\ 

Taking Vh = Wh above and using the Poincare inequality, we know 

\wh\i < -p°||o- 
Also, by the standard error estimate, we have 

Now, taking v = — in (|45|) and using 1)441) and the duality argument, we obtain 

Wp'-pUI = {Vw,V{p'-pI)) 

= {V{w - Wh),Vip'' - pD) + {Vwh, V(/ - pD) 
= {V{w - Wh),V{p^ - pI)) + {V{lihs - s), VK - w)) 
+{V{Iih~s-~s),Vw) 

< \w - Wh\l |p° -Pl\l + \^hS - s\l \wh - W\i + \IlhS - s\l-ao\w\l+ao 

< h^-0\\p^-pl\\o + /i^+°°|s|2||p° -P°J|0 , 

which leads to the desired L^-norm error estimate: 

WPs - P% < - pUo + \\s - lihSWo < + h^\~s\2 < h"^^. ^ 



Remark 5.1 Following the proof given in one can improve the results of the previ- 
ous Lemma. Indeed, one can derive the estimates \ps—p^\i ^ cLnd \\ps—p^\\o ^ 
with slightly more restrictive assumptions on the mesh. 



5.2 Approximation of the singular part 0^ 

In order to approximate the singular part (j)s in the decomposition u^ = + , we 
recall (cf. [0|) that 4>s G Hq{uj) solves the elliptic problem (|22|) and has the following 
decomposition: 

0, = (^ + /3>p, 4>^H''{lo), /3* = i||p,||2, </,^ = sin(a(A) . (46) 

TT 

Using ()22() . we see that (j), satisfying (f) = —(3*(j)p on duj, solves the variational problem: 

{V^,Vv) = {ps,v) yveH^iio). (47) 
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The next step is to consider the finite element approximation of (p in V^: 

4>h = -PlT^h(t>p + (Pt 

where tt/j is defined as in (|in|l . is computed using /3/^ = — / {Ps)'^d^^, and (f>1 G Vq 
is the solution to the problem: 

iV4>h,Vvh) = {p':,Vh) yvheVo\ (48) 
Then we propose to compute the finite element approximation of (pg by 

0^ = 4 + pi(Pp . 

Below, we derive the error estimates for this approximation. 
Lemma 5.2 The following error estimates hold 

\^s-<P%<h, \\4>s-(p%<VJ^h. 
Proof. We first estimate the error cp — cp^. Subtracting ()48() from H47|) yields 

(V(0 - 4>h),Vvh) = ips-P^vh) yvhe < , 

thus we obtain for any Wh G satisfying Wh — (ph ^ Vq^ 

- Wh\l = \(p- Cphll + \(t>h - Wh\l + 2{ps - p^, (ph - Wh), 

which with the Young inequality and the Poincare inequality gives 

\^-4>h\l < \4>-Wh\l + 2Cp\\ps-p%i\4>-k\i + \4>-Wh\i) 

< 2\4>-Wh\l + \\4>-Ml + \\Ps-Ps\\l (49) 

Noting that (p = —P*(pp on div, so (5lllh(f> = l3*(ph on duj. Let Wh = (3^11^4)/ /3* , then 
Wh — 4>h ^ y^o - With this Wh, we derive from 1)491) and Lemma l5.ll that 

\4>-k\l < h^ + m-^\f3*4>-(3h^h^\l 

< h' + \p*-f3i\'\4>\l + \/3l\'\4>-iih4>\l (50) 

But using the definitions of f3* and /?^, we have 

\(3^-Ph\ = -\\Ps\\l-M <\\Ps-p%<h'-'. (51) 
vr 

It follows from ()5()j) and the property (p G ff^(a;) that 

\4> - 4>h\i < h. 

This with ()51() and the decompositions of </)s and 0^ gives the desired iJ^-norm estimate: 

\cps-(t>s\i < \4>-4>h\i + \P'' -Pl\\4>p\i<h. 

Finally, by noting that both (pg and <pg vanish on 71 and 72, we can apply the Poincare 
inequality to the function (ps — (p^ to get 

\\<Ps-(p%<C'p\<Ps-(P%. 

Then the desired estimate on ^4>s — 4^s 1 1 (I follows from 

Us - cp^Al = \cPs - 4>'l\l + 1^ Us - 4>X <h^ + l^h\ o 
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5.3 Approximation of and in decomposition 

Noting that and solve the coupled system and (|32|) , it is natural to formulate 
their finite element approximations as follows: 
Find u^^ G Vq and c'^ £ M} such that 

a^{u%v) + cla^{<Ptv) = {f.v) Vt; G (52) 
{\\pX + + H {u%p^) = {f,p^) , (53) 

where (j)^ and are the finite element approximations of (ps and ps, see Subsect. l5.1l5.2l 
However, this formulation requires solving a coupled system, and it poses some 
difficulty in getting the error estimates as it does not fall into any existing saddle- 
point-like framework. Instead, we are going to propose a more efficient approximation 
which enables us to find G Vq and separately. In fact, we can use the formula 
to first find c|^, and then use to find tt|j^ G Vq. This leads to the following 
algorithm to find G Vq and c|j. Let C* > be a fixed constant. 

SCM Algorithm for finding G Vq and c'^ G 



Step 1. Find G Vq such that 



a^{z';,v) = {f,v) ^veV^. (54) 

Compute as follows: 



WPsWo 

and 

Ji 



-'= ^^ .ri/'^ if VJ^<C*h-^o. (55) 



= if ^>C*h 2-°o . (56) 
Step 2. Fmd G Vq* suc/i i/iai 

a^{u''^,v)+c''^a^{(t>'',,v) = {f,v) Vt; G < . (57) 



Remark 5.2 In practice (see I15f ). the conditions il55H56)] mean that only a few coef- 
ficients (c^)^ neec? to be computed, with respect to the total number of Fourier modes. 

Below, we shall derive the error estimates on (c^ — c^) and {u^ — u'^). Recall the formula 
dSni) for c^: 

_ {f-fiz^,ps) 

WPsU ' ^ ^ 

where z^ = A~^f G -f^o('^) solves 

a^{z^,v) = if,v) yveH^iu). (59) 

Clearly Zf^ = Uf^, the solution to the equation (PU)) . But a different notation z^ is used 
here for convenience, since the numerical approximation z^ is derived with the standard 
piecewise linear FEM. 
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Lemma 5.3 For the solution to the problem h5y\) and its piecewise linear finite 
element approximation in \54^ , we have the following error estimates 

Ik^-^llo < /^"l/llo, (60) 



hn^ < /i2-o^-o-i(i + ^/,)2||^||^^ (61) 



while for the coefficient c^ in i58\) and its approximation in k55\} . we have 

|cm - clJI < + ^hf + h) ll/llo . (62) 

Proof It follows from and (EH) that 

CLfiiz/j, ~ Zfxi Z/j, — z^) = a^(z^, z^ — z^) = (/, z^ — z^). 

This implies 

I^M - 4\i + ^' II^M - 4\\o < ll/llo \\Zf, - ^^llo, 

thus (|6fl|) follows by the Young inequality. 

We next show Again it follows from (|54jl and (|59|) that 

But, by standard interpolation theory, we know that 

\z^-IlhZfj,\i < h"°\z^\i+ao, and ||z^ - XI/jZ^Ho < l^^li+ao- 
Therefore, we reach 



z„ 



),-z^Ja<{l + VJ^h)h''°\z^\i+^,. (63) 
On the other hand, for any g G L'^{uj), define w S ^^o(<^) such that 

af,{w,v) = {g,v) yveH^iuj). (64) 
Using the duality and (|64|) . we have 

II hu (^M-^M'5) a^(w^,2;^ - ^/J) 

lF/i-^ul|0 = sup rr-rf^ = SUp rr-r 

geL^ioj) \\9\\o geL^iuj) \\9\\o 

a^{w -UhW,Zf, - Z^) \\W -UhW\\a\\Zf, - z'^Wa 

= sup — — < sup — — . 

c/eL2(c^) ll/llo 5eL2(a;) llsllo 

Using the interpolation result and the same derivation as in and the a priori 
estimate ()26p (with u and / replaced by w and g), we obtain 

II hii ^ /l^°°(l + V/^/l)^kll+aol^A«|l+ao 

IFm-^uIIo ^ sup ^ — 

geL2(<^) ||5||0 
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which proves H61|) . 

It remains to prove ()62|). We have from H55|) and H58() that 

wpsU mi 

For /i, we have from Lemma lS.ll that 

\h\<h\\f\\o. 

For I2, we further write it as fohows 

Then using estimate (|61|) and Lemma l5.ll we can derive 

\l2\ < /X 11^^ - z;:||o + h ll/llo < (/.2-o^«o(i + ^;,)2 + f^-^ ii^ii^ _ 

This with the estimate of Ii gives 1)62^ . o 

In the rest of this Section, we shah estimate the error between the solution to the 
ehiptic problem ^lU^ and its SCM approximation n|j. We note that the decomposition 
of is equal to: 

= + c^(j)s = Uf^ + c^{(f> + (3*(j)p) . (65) 
So, we propose its SCM approximation of the form: 

< = ^^1^ + ^''s=u''^ + (4 + l3l<Pp). (66) 

We shall derive the error estimate on — u^- Let us start with the estimate of 
{Ufj, — u^)- We have 

Lemma 5.4 The following error estimate holds 

\\n,-uX<^{k'\\f\\l + \c,-c'lf). 
Proof. Subtracting (|52|) from ^?>1\ we have 

a^{u^ - H'l^Vh) + c^a^{(j)s,Vh) - c^^a^{4fl,Vh) = M Vh G Vq ■ 
Using this we obtain for any Wh G Vq, 

\\u^, - WhWl = - u^Xa + H^J ~ '^^Wl + 2cja^('/'s > - Wh) - 2c^a^(0^, - Wh), 
which implies 

< II^M-^/^lla + 2c;.a;,(</., -</.^,n|i-u;/,) + 2(c^-c|j)a;,(</>^,'i2'^-u;/,) (67) 



< 



h\\ llr.h II I o l„ „h\ \\j.h\\ \\~,h 



\Uf^ - WhWi + 2 \C^,\ Us - Clla Wu"^ - WhWa + 2 - C^^j Wja " Wh\\a ■ 
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Now, there holds ||(/'s||a- \\4>s - (t>s\\a < lltAslla < ll'/'slU + ||<^s - '/'^lU- Using LemmaEl 
and ||i;^<i||a = l^sli+A* ll'/'sllo' '^^ ll<^s lU ~ \//^ Il0s||o- Using the interpolation results, 
we obtain 

\\u^ - IlhUfj.\\l < \Ufj. - UhU^ll + fj, \\u^ - II/iM^llo < \u^\l, 
thus letting Wh = H/jU^ in 1)6 7p and using Lemma 14.21 we derive 

'\U„. - u'l.Wl < h^\Ut,\l + ^/J^h'^\\f\\o\Uf^\2 + ^/J^h\C|,-c'|^\\Uf,\2 



< VJ^ih^\\f\\l + K-c',\'). o 

Theorem 5.1 Let be the solution to the equation \2U\) and be its finite element 
approximation given in if) 6}) . Then the following error estimate holds: 

3C > such that V/u, ||n^ - u'^\\a < C fxh\\f\\Q. 

Proof. It follows from H65|) and ()66() that 

Then we obtain, using Lemmas 15.41 15.21 and 14.21 that 

\K-uX < 3{||S^-n;:||2 + |c^|2||,/,,-0,^||2 + ||</,^^||2|c^_c;^|2| 

< ^h'\\f\\l + ^,\c^-c^^\\ 

To prove the desired estimate, we need simply 

\c,-c';,\'<f^h'\\f\\l (68) 

1 

First consider the case (|5()|) . i.e., ^/JI > C* h ^-^o . This condition is equivalent to 
Then (|68|) comes directly from this condition, = and (|25|) as follows: 

A|2 _ „2 ^ ,.oo-l|| -p||2 ^ ,. 1,2 /i,-2 ,,ao-2Ml ^||2 ^ ,. t 2 ii /^||2 



Faj-C^I — Cf, rZ H- WJWOrZH'"' H- )\\J WO rZ H- WJWO- 

1 

For the remaining case (|55|) . we have -y//! < C* h ^~°o , or /i^ < /i^'^^""'. On the 



one hand, since ao < ~ /i^""o < 1. On the other hand, since 2ao — 1 > 0, 

^4ao-2 < ^-(2ao-i)(2-ao)_ g^j; Qj^g infers from and these inequalities that 

To conclude, (|68jl follows from this and the fact that, as ao s]^) ![> the exponent of 
is bounded by 

2ao - (2ao - 1)(2 - ao) = 2a^ - 3ao + 2 = 1 + (2ao - l)(ao - 1) < 1 . o 
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6 Fourier Singular Complement Methods 



In order to define the numerical part of the Fourier Singular Complement Method, let 
us prove a result which can be viewed as the mathematical foundation of the FSCM, 
from the Fourier point of view. It allows to recover (|4l6jl . for sufficiently smooth right- 
hand sides. 

Let u be the solution to the Poisson problem and Uk be its Fourier coefficients 
in (|T6|) . By LemmaEIH we know that Uk{xi,X2) solves the 2D problem (|17II18|) . And 
using (pS)) we can decompose Uk as follows: 

Uk = Uk + Ck 4>s (69) 
where Uk G H'^i^^) n H^{uj) and (/j^ G H^{^) solves (|^ . 

Lemma 6.1 Let f G h^{Q) n h^{Q), and u G Hq{Q) be the solution to (CJ). Then 

u = u + -f{x3)(l)s, withuG H^{n)nH^{n), -f e H^{]0,L[)nH^{]0,L[). (70) 

Proof. Let {Uk)k be the Fourier sequence of u. Recall that (Uk)k converges to u in 
Hq{Q,), and {AUk)k converges to — / in L'^{Q). From (|69|) . let us split the Fourier 
sequence into regular and singular parts, as 

K , K , 

f^x = f^x + 7x(a;3) <As, with ?7i^ = 2^ sm — xa, -jKyxz) = ^^Cfcsm— X3. 

fc=l k=l 

We shall prove below that {'~^k)k converges in -?/^(]0, L[) n i^o (]0, L[), and {Uk)k con- 
verges in H^{^)r\Hl{^). 

As far as the singular part is concerned, from ()14() and the bound on \ck\ in 

00 

Lemma 14.21 we obtain that /c'^|cfc|^ < 00. Since we are dealing with the ID Fourier 

k=l 

sequence {'^k)k (with sine functions), it is well-known that it converges to a limit, sub- 
sequently called 7, in -ff^(]0, L[) n -f^o(]0, L[). Then, one finds that {'^k 4's)k converges 
to 7 03 in Hq{Q), and that {A{'yK 4's))k converges in L^(r2), to 7" (/>s — 7Ps. 

For the regular part, we note that since there holds Uk = Uk — 1k4>s-, {Uk)k 
converges in ifQ(O), to a limit called u, which is equal to 

■U = M - 7(^s. 

Moreover, {A.Uk)k converges in L^(r2), to Au. 

To conclude the proof, one has to establish that u is an element of i/^(J7). From 
Corollary ESI we know already that d^u is in H^{^). So one has to check that dijU is 
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in L^(0), for i,j G {1,2}. But this follows from the estimate on |nfc|2 in Lemma 14.21 
and on the expression of the second order partial derivatives of Uk , that is 



k=l 



Remark 6.1 In the more general case, i.e., f G L^(S7), one gets only a convergence 
of {ik)k H^~'^(\^-,L\), see 11 1^ . This precludes a convergence of the singular part 
in the desired Soholev spaces, i.e., H^{Q) with L^(Q) Laplacian. 

In order to build the numerical schemes which completely define the FSCM, we intro- 
duce u^{xi,X2) the SCM approximation to Uk{xi,X2). It is the same as in H66|). but 
with /i replaced by /c^vr^/L^, that is, 

We then rephrase the 2D SCM Algorithm H54I57|) . This gives 
Step 1. Find G Vq such that 

ak{ztv) = {f,v) yveVo^ (71) 

Compute as follows: 

4 = hr2 if k<C*-h-^; (72) 

l|P.llo ^ 

and 

4 = if k>C*-h~^. (73) 

TT 

Step 2. Find G Vq such that 

ak{vll,v)+clak{(t}'l,v) = {f,v) "iv^V^. (74) 
As mentioned already, only a few coefficients {c^)k are actually computed. 

Following ()19|) . we finally define the FSCM approximation to the solution n to 
as follows: 

N 

U'}:i{xi,X2,X^) = ^U^(xi,X2)sin— X3 . 
k=l 

Then we have the final error estimate below 

Theorem 6.1 Assume that f G hl{Vl) n h'^{VL). 
The following error estimate holds: 

||V(n - UU\l2(^^) <{h + N-'){\\f\yin) + Wdssfh^^n)} ■ 
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Proof. Using the Fourier expansion of u and the definition of C/^, we have, cf. HlOj) . 

k=l 

4E(iiv-.iig + (x)'ii-^iio 

k>N 

=■■ I1 + I2. 

According to Lemma 13.41 we derive 



h 



k>N 

L . o \ T ,0/11 1 1 o - AjTt , 



2|L. I|2 



For Ii, we have 



fe>Af 



^ h\\2 



2 

k=l 



According to Theorem 15 . 1 1 we have 

\\uk-4\\l<k'h'mi. 

Using this and ()14() . we obtain the estimate of Ii: 

N 

h<h'Y.^'\\fk\\l<h'\\d,,f\\l,^^^, 

k=l 

which, together with the previous estimate of I2, leads to the desired error estimate, o 

7 Conclusion 

The optimal convergence rate of the FSCM in prismatic domains, has been proven for 
the Poisson problem with homogeneous Dirichlet boundary conditions. Assuming that 
the right-hand side / is slightly more regular than / G L^(r2), i.e., that / belongs to 
n hl{Q), the convergence rate of the FSCM in if^-norm is like 

\\u-U^\\i <Cf{h + N-^), 

where h is the 2D mesh size, and is the number of Fourier modes used. 
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The same result also holds for the discretization of the Poisson problem with a 
homogeneous Neumann boundary condition, or for the Poisson problem with non- 
homogeneous boundary conditions, provided there exist sufficiently smooth liftings. 

Further, it is no difficulty to consider the case of a prismatic domain with several 
reentrant edges, i.e., oj with several reentrant corners. 

As far as the assumptions on the right-hand side / are concerned, a few remarks can 
be made. It seems that, in a prismatic domain Q. such as the one we considered here, 
the boundary condition on the bases was omitted in [2]. Nevertheless, this condition 
does not exist in the case of an axisymmetric domain, see jl4j . nor in the case of an 
infinite cylinder. In other words, / € h?'{Vl) is enough in those types of domains. In 
the case of a Poisson problem with Neumann boundary conditions, one has to replace 
the vanishing trace conditions at the bases by the familiar d^f = at the same bases. 

As mentioned already, this paper is the first part of a three-part article |141ll5j . In 
the companion paper the FSCM is analysed theoretically and its numerical ap- 
proximation is built, in axisymmetric domains with conical vertices and reentrant edges. 
There are two difficulties which are inherent in this class of domains. The first one is 
the weights, which have to be introduced in the 2D sections. The second one is the 
addition of sharp vertex singularities, which have to be taken into account separately. 
In jl5j . the FSCM is analyzed from a numerical point of view (complexity, implemen- 
tation issues, numerical experiments, etc.), and it is compared to other methods, such 
as mesh refinement techniques, or variants of the FSCM (2D SCM with the A-approach 
|13j : 3D discretization of the regular part, etc.) in prismatic or axisymmetric domains. 
In particular, the use of the FFT to aproximate the sine functions in 2:3 is motivated 
and justified there. 

As noted in Remark 14.11 one can apply the same theoretical and numerical tech- 
niques to the 2D heat or wave equations, with any L^-smooth (in space) right-hand 
side. For these PDEs, the singular functions ps and (ps do not depend on the time- 
step. 

Finally, the results, can also be viewed as the first effort towards the discretization 
of electromagnetic fields in prismatic domains, with continuous numerical approxima- 
tions, the importance of which is well-known, cf. [5]. As a matter of fact, the SCM 
developed in [SI [71 for 2D electromagnetic computations can be generalized, based 
on the results obtained here. 
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